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Abstract 

Time evolution is formulated and discussed in the framework of Schroder's functional equa- 
tion. The proposed method yields smooth, continuous dynamics without the prior need for local 
propagation equations. 
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I. INTRODUCTION 



From a lattice of time points, as is usually discussed in dynamical systems, is it possible 
to obtain continuous, or even smooth, time evolution without ambiguity? Under certain 
circumstances, the answer is yes, through a holographic interpolation process involving func- 
tional methods. 

In this paper, we utilize and illustrate Schroder's functional techniques [1] to produce 
continuous trajectories out of discrete nonlocal recursion laws (orbits) for less self-evident 
dynamical systems, with fixed points at x = 0. In particular, we build exact continuous 
iterations of 2x (1 + x), and approximations to continuous iterations of xexpx, and their 
inverses, and we discuss the notable features of the resulting evolution. Typically, smooth, 
analytic interpolates of the discrete recursion rules that act as the boundary data for the 
process (hence our use of the term holographic) are obtained. These methods are of possible 
use to anyone pursuing practical applications of functional equations 0, and are also likely 
to be relevant to more recent theoretical developments 0]. 



II. THE METHOD 

Consider an evolution trajectory x (t) of a 1-dim system, e.g. specified by a local, time- 
translation- invariant law (cf. energy conservation). 

dx(t) /alt = v(x(t)) . (1) 

One may integrate this to obtain the trajectory as a family of functions of the initial data, 
indexed by the time, 

x(t) = f t (x{0)) . (2) 
For any given time, we may scale t to consider a unit time increment At = 1, so that 

x(l) = A(x(0)) . (3) 

But then, for time-translation-invariant systems, further evolution obeys 

x(t + l) = f t+1 (x(0)) = f 1 (x(t)) , (4) 

i.e. x (t + 1) is the same function of x (t) as x (1) is of x (0). For notational convenience in 
the following, we will often use x (0) = x. 

Now suppose, for reasons dictated by physics applications, that the time-local evolution 
law ([1]) is not specified, but an explicit, nonlocal, discrete propagation function f\ is avail- 
able, as in 03]). We pose the question: How does one obtain the complete, continuous 
trajectory x (t) = f t (x) without benefit of the local relation? 

Of course in principle, it is straightforward to compute iterates of fl3]) on an integral lattice 
of time points, t = ■ ■ ■ , —2, —1, 0, 1, 2, 3, • • • , to obtain "the splinter of x" (for example, see 
@) , i.e. 

x(2) = fi (fi (x)) = f 2 (*) , 

x (") = fi (/i • • • (fi (*))) = fn (x) , (5) 
x(-l) = /f 1 (x) = U (x) , 
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etc., assuming the domains for the various functions overlap properly. Thus, x = 
/-l (/i (a?)) = or more generally, x(k + n) = fk(fn(x)) = f n (fk(x)), asso- 

ciative and commutative composition. From this lattice, upon selecting derivatives at the 
lattice points, an infinity of interpolating functions could be produced, say graphically, to 
obtain a continuous trajectory. This approach permits easy visualization of the trajecto- 
ries, but is ambiguous, and does not place proper emphasis on analytic properties of the 
solutions, so it will not be pursued. 

Instead, given fx (x) we will use here the theory pioneered by Ernst Schroder [l| to 
construct an analytic f t (x) around a fixed point of fx (x). Without loss of generality, we 
take the fixed point to be x = 0. 

Schroder's construction of ft (x) amounts to building all iterates of fx (x), including frac- 
tional, negative, and infinitesimal t, based on his eponymous functional conjugacy equation 
involving the auxiliary function \I> Q], 

s*(x) = *(f 1 {x)) , (6) 

for some constant s / 1. With the origin a fixed point of fx, i.e. fx (0) = 0, it follows 
that \& (0) = 0, and if ty' (0) ^ 0, then s = f[ (0). The inverse function satisfies Poincare's 
equation, 

vp-i ( sx ) = f x ( x )) . (7) 
Upon iteration of the functional equation, \1/ acts upon the splinter of x to give 

S n * (X) = V {fn (*)) = * (fl (/I""" (fl (*)))) • (8) 

Now, the point to be stressed is that this formula naturally yields a continuous interpolation 
for all non-integer n, 

s t *(x) = *(f t (x)) . (9) 

So, to produce the full, continuous trajectory, we solve for Schroder's auxiliary function 
*&(x), and construct the inverse function Having done so, this yields x (t) as a 

functional similarity transform of the s* multiplicative map. From Q, 

x (t) = f t (x) = ^ l (s* * (x)) . (10) 

In a suitable domain, this gives the general iterate for any t, analytic around the fixed point 
x = 0. 

Moreover, this solution manifestly satisfies the requisite associative and abelian composi- 
tion properties for all iterates and inverse iterates. That is to say, f^+tn ( x ) = fti [ft 2 ( x )), 
hence x (tx +t 2 ) = ft t (x (£2)); as required for time-translationally invariant systems. Some 
specific cases are as follows. 

f 2 (x) = V- 1 (s 2 * ( x )) = ( a * ( fl ( x ))) = fl (/l ( x )) , 

fx (X) = ^ l (S 1 * ( x )) = «r-l ( S V2 y ( s l/2^ = /i/2 ( /i/2 ) (n) 

f (x)=X = {S- 1 * (^ (s 1 * (X)))) = /_! (/1 (z)) , 

etc. However, it is crucial to note that in the limit s — > 1, all iterates and inverse iterates 
lose their distinction and degenerate to the identity map, f (x) = x, and the method fails 
as written. For this reason, if f[ (0) = 1, one augments fx (x) in Schroder's equation to 
sfx (x), and takes the s — > 1 limit only at the very end of the calculations, if it makes sense 
to do so. 
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III. SCHRODER'S EXAMPLE 



For a very elementary illustration of the technique, consider Schroder's early example of 
a recursive evolution law, 

f 1 (x)=2x(l + x) , (12) 
so fi (0) = 0, and s = f[ (0) = 2. Schroder's equation §6§ is then solved by 

# (x) = \ In (1 + 2x) , (x) = \ (e 2x - l) . (13) 



■rih LU) ~((n 2,-r - 1) , (ID 



This results in 

X (t ) = f, (x) = 

which indeed obeys f h (f h (x)) = f tl+ t 2 {%)■ In particular, f x (x) = |((1 + 2a;) 2 - l) = 

fi/2 (/i/2 (x)), f 1/2 (x) = \ ((1 + 2x) V ~ 2 - l), /r 1 (x) = /_! (x) = | (v/T+2^ - 1), etc. 

We may visually appreciate the evolution described by (fl4"j) through a surface plot, where 
the left-right axis is the initial position, x, the depth axis is the time, t, and the vertical axis 
is x (t). The surface flows from the identity map at t = to the given discrete propagation 
function x i — )■ /i (x) at t = 1. The origin is invariant since it is a fixed point of all the 
functions /t(x), for every value of the t index, as is the point x = — 1/2 in this simple 
example. Also note, here we require the initial x > — 1/2 to have real x (t) for all t. 




Figure 1: x (t) plotted versus t and initial x for Schroder's example. 

The full sweep of the surface describes not just one trajectory, starting at a single value of x 
and moving into the page with t, but rather it encodes the evolution of all trajectories from 
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the initial domain x > —1/2. The construction is holographic in the sense that the interior 
of the surface is completely determined by the behavior on the time boundaries, namely, the 
identity map at the front edge and f\ (x) = 2x (1 + x) at the back edge. An animation of 
sequential, fixed time, vertical slices through the surface would perhaps show more clearly 
how the initial straight line of data evolves to the final fi [Til ] . 
The velocity profile following from ( |14p is 



v(x(t)) 



dft 



x 



dt 



[l + 2x (t))ln(l + 2x (0)lnV2 . 



(15) 




Figure 2: Velocity profile for Schroder's example. 

Had this velocity been available ab initio, it could have been integrated in the usual way, 
i.e. t = ^y, to obtain the trajectory (1T4"]) moving into the page for each x as shown 

in Figure 1. But here the starting point was the discrete time step ( TT2l) . with v (x (t)) an 
emergent feature flowing from the formalism 0]. 

Moreover, assuming the motion is governed by a Lagrangian, with L = |?m> 2 — V(x), 
and that we are dealing with this system at a fixed energy, the result for the velocity profile 
immediately gives the ^-dependence of the corresponding potential. Namely, 



V(x) 



-mv (x) + constant 



(16) 



In this sense the functional method of determining the evolution surface is a technique to 



solve an inverse problem [13l.ll4j: Given the fi (x) "scattering data" for a finite time step (as 
opposed to the usual infinite time step in an idealized scattering process), we may determine 
an underlying V (x) upon making certain analyticity assumptions about the solution. For 
the Schroder example, the essential features are contained in a plot of the effective potential 



-v IX 



((1 + 2x) ln(l + 2x) lnv^)' 
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V(x) 




Figure 3: Potential for Schroder's example. 



Note this V is unbounded below, and more importantly, the aforementioned fixed points 
at x — —1/2 and x = are simply unstable stationary points of the potential. We submit 
that this is a common feature of the functional method, when the results are expressed as a 
potential function. 

In general, a formal solution of Schroder's equation is predicated on the existence of 
infinite iteration limits of fi (x), as discovered by Koenigs [7|. (The technique itself is 
perhaps more familiar in the context of the Poincare equation when written as a nonlinear 
finite difference equation.) When the splinter is "approximately geometric" (for example, 
again see [5|) we have 

(x) oc lim ^M. . (17) 

7V->oo s N 

In the above example due to Schroder, fTH|) . we actually require the limit as t goes to negative 
infinity, 

9(x) = ±\n(l+2x)= lim (18) 

Z N— >oo Z 

as it is the splinter of /_i that is approximately geometric. 



IV. ALL ITERATES OF THE RICKER MAP 

We now consider the iterates of the Ricker map [§], 

fi (x) = xexp (x) , (19) 
or equivalently, of its inverse, the Lambert function [9J, 

/_i (x) = LambertW (x) = ^ x n . (20) 

n=l U ' 

The latter series is convergent for \x\ < 1/e. Both f\ (x) and /_i (x) have a unique real 
fixed point for x = 0, so Schroder's conjugacy equation can be solved for the auxiliary ^ 
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with the same fixed point, in power series expansion around that point. Since f[ (0) = 1, 
as explained above, we will modify (|T9|) and (I20|) to 

/i (x) = sxexp (x) , /_i (x) = LambertW (x/s) , (21) 

thereby introducing a second fixed point in fi at x = — In s, but allowing for a possible limit 
s — )■ 1 at the end of various calculations, if sensible. 

The resulting Schroder's equation looks deceptively simple, 

s^j ( x ) = * (sxe x ) . (22) 

Here it is implicit that ^ really depends on two variables, both s and x, so (122]) is actually 
s\I/ (x, s) = \I/ (sxe x , s). In the following, this additional dependence on s will be understood 
but usually not displayed. It follows from (|22|) that s^/' (x) = se x (1 + x) (sxe x ), along 
with all higher derivatives with respect to x given by Faa di Bruno's general formula, hence 
an explicit series solution for \l/ about x = is straightforward to construct, in principle. 

The corresponding Poincare equation for the inverse function, $ = \J r-1 , while nonlinear, 
is also relatively simple in appearance (lot ]. 



- $ (sx) = $ (x) exp ($ (ar)) . (23) 

Again, it is implicit that $ depends on both s and x. 

With the normalization choice (0) = 1, the auxiliary function and its inverse are given 
explicitly to O (x 5 ) by 

, > 1 2 1 3s + 1 , 1 16s 3 + 8s 2 + lis + 1 . 

w It I = t — t -I- t — nr. 

K) (s-l) 2(s-l)(s 2 -l) 6 (s-l)(s 2 -l)(s3-l) X 

1 125s 6 + 75s 5 + 145s 4 + 146s 3 + 53s 2 + 31s + 1 5 , 6 , 
+ 24 (s - 1) (s 2 - 1) (s 3 - 1) (s 4 - 1) X + °y x ) ' ^ 



! 1 2 1 3 + s 3 1 16+ lis + 8s 2 + s 3 4 

[X) ~ X + (s - i f + 2 (s - 1) (s 2 - 1) X + 6 (s - 1) (s 2 - 1) (s3 - 1) X 

1 125 + 131s + 145s 2 + 106s 3 + 53s 4 + 15s 5 + s 6 , ^ , fiN 

H ; ; x + O (x b ) . (25) 

24 (s - 1) (s 2 - 1) (s 3 - 1) (s 4 - 1) V ; V ; 

For purposes of illustration, we plot O (x 10 ) approximations for the auxiliary function and 
its inverse versus x for various values of s, to obtain the surfaces below. We wish to convey 
only qualitative behavior at this point, not detailed structure, and to point out that for 
values of s > 1 the auxiliaries have many features similar to those for Schroder's example, 
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1.5 

Figure 4: \1/ contour surface to O (x 10 ) for Figure 5: \E' _1 contour surface to O (x 10 ) for 

-0.5 < x < 0.5 and 2.0 > s > 1.5. -0.5 < x < 0.5 and 2.0 > s > 1.5. 

We stress that these are only approximate representations for the auxiliaries, based on the 
explicit series to O (a; 10 ) . We do not give the O (x 10 ) series explicitly for all s, although 
they are straightforward to obtain. For example, at the back and front edges of the plotted 
surfaces, the series are given numerically by 

* 0) | s=2 = x - l.Ox 2 + 1. 166 7a; 3 - 1. 452 4x 4 + 1. 873 4a; 5 - 2. 470 8a; 6 + 3. 308 5a; 7 

-4. 478 8a; 8 + 6. 113 3a; 9 - 8. 398a; 10 + O (a; 11 ) , (26) 

( x ) | s=2 = x + x 2 + 0.833 33a; 3 + 0.619 05a; 4 + 0.424 21a; 5 + 0.273 6a; 6 + 0.168 26a; 7 

+ 9. 952 9 x 10~V + 5. 698 9 x 10" 2 x 9 + 3. 173 4 x 10~ 2 a; 10 + O (a; 11 ) , (27) 

* ( X )L=3/2 = X ~ 2 - 0x2 + 4 ' - 10 - 049x4 + 23 - 402x5 - 55 - 143x6 + 130 - 95x7 

- 312. 72a; 8 + 749. 86a; 9 - 1803. 8x 10 + O (x 11 ) , (28) 

^ ( x ) | s=3/2 = x + 2.0x 2 + 3. 6a; 3 + 6. 049 la; 4 + 9. 667 3x 5 + 14. 861a; 6 + 22. 142x 7 

+ 32. 145x 8 + 45. 656a; 9 + 63. 633x 10 + O (x 11 ) . (29) 

Also, we have avoided s = 1 in the plots of Figures 4 and 5 since the explicit series results 
exhibit an expected singular behavior as s — > 1 for each of \1/ and \J r-1 , considered sepa- 
rately. This behavior is foreshadowed by the growing coefficients in the numerical series for 
tt(ar)|, =3/2 and ^ (x)| 5=3/2 . 
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However, when \l/ and \I/ 1 are composed as in (fTUj) . the result for /( (x) is well-behaved 
even in the limit s — > 1. Take s = e e , and expand in powers of e, to find 

/t(*)|«.= *- 1 («'*(ar))| -==e . (30) 
= (l + te + 0(£ 2 ))x 

+ U + ^ (-1 + 3t) te + O (e 2 )J a; 2 

+ Q (-1 + 2t) t + ^ (-1 + 2t) 2 te + O (e 2 )^j x 3 

+ ( — (5 - 15t + 12t 2 ) £ + — (-7 + 35t - 56t 2 + 30t 3 ) te + O (e 2 ) ) a; 4 

+ f ^ (-2 + 3t) (5 - 12t + 8t 2 ) t + (50 - 315t + 673£ 2 - 621t 3 + 216t 4 ) te + O (e 2 )^j x 5 
+ 0(x 6 ) . 

Despite what one might naively expect from the form of fi in (j2Tl) . the dependence of /t 
on s is certainly not multiplicative, a point already borne out by /_j in ff2TT) . Of course, 
as i — >• all e dependence (i.e. all orders in e) must disappear to yield the identity map 
fo (x) = x, and indeed (13171) reduces accordingly. 

This last result permits us to compute, at least to O (x 5 ), the initial velocity profile in 
the limit s — > 1, to find 



A better numerical approximation keeps all terms to O (x 10 ), to give 



(31) 



v (?) \ s=1 = x 2 - 0.5a; 3 + 0.416 67a; 4 - 0.416 67a; 5 + 0.445 83a; 6 - 0.480 56a; 7 

+ 0.501 12a; 8 - 0.491 63a; 9 + 0.452 15a; 10 + (a; 11 ) , (32) 

as plotted here. 

v(x) 



Figure 6: Initial velocity profile to O (x ) as s — > 1. 
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Moreover, from time-translational invariance, we automatically obtain power series approx- 
imations to v (x (t))\ s=1 for all t just by substitution of x (t) in either of (|3"T|) or (I321) . We 
may also visualize the effective potential for this problem from a plot of — v 2 (x), as discussed 
previously in the context of Schroder's example. 



V(x) 



Figure 7: Effective potential to O (x 10 ) as s — > 1. 

As in that previous example, we see the fixed point at x = is a point of unstable equilibrium. 

A similar calculation produces the initial velocity profile for other s ^ 1. Again to O (x 5 ), 
we find 



v (x, s) = lim = li m A^- 1 (V ^ ( x s ) s ) 



(In s) I x + 



1 3 1 
^ + - 



3s + 2 



^s 2 + 4s + 3 



j 2 -! 



2 s 2 -l s 2 + s + l 



-X 



3 (s 4 - 1) (s 2 + s + 1 



-X 



(33) 

> 1 O (x 6 ) 



Here we have made explicit the dependence on both x and s. Of course in the limit s — > 1, 
f (x, s) remains finite and reduces to the previous (13"T|) . We plot the initial velocity surface 
versus x and s. 



v(x,s) 




1 1.0 



Figure 8: Initial velocity contour surface, to O (x 5 ), for — 1 < x < 1 and 1 < s < 2. 



10 



As before, from time-translational invariance this function automatically also yields 
v(x(t),s) for all t just by substitution of x(t) in the power series ( 133|) . In addition, 
the potential surface is effectively given by V (x, s) = —\m (v (x, s)) 2 . 



V(x,s) 




1 1.0 



Figure 9: Effective potential surface, to O (x 5 ), for —1 < x < 1 and 1 < s < 2. 

For fixed s slices, this potential surface again shows unstable equilibria at the fixed points, 
x = and x = — In s, as was the case for Schroder's elementary example. 

To visualize the evolution, we again plot x (t) versus t and initial x, as we did for the 
closed form results of Schroder's simple example. We find the same general features. 
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0.0 _! 



Figure 10: x (t) = lim f t (x), to O (x 10 ), plotted versus t and initial x. 

s— >l 
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We only display the limit as s — > 1. Other s give similar surface geometries. As before, 
the front edge of the surface is the identity map, and the back edge is the discrete propa- 
gation function f\. Here f\ (x) = xexpx. Continuations out of the page would give the 
corresponding inverse functions, /_ t , including /_ x (x) = LambertW (x). We reiterate that 
this evolution surface was not obtained by the standard method of integrating velocities for 
individual initial x to produce the lines that go into the t depth of the mesh, but rather by 
the use of functional methods to produce the continuous, inward flow of complete left to 



right "time slices" of the surface [11 



With these qualitative images in view, we turn to discuss in more detail the series solutions 
for the auxiliaries and the composite functions ft (x) that follow from f\(x) = sxexpx. 
Generally, the series for \l/ and are of the form 



71=1 V ' \k=l / 71=1 V ' \fc=l 

The Pn are polynomials in s, and can be obtained explicitly by recursion from pi — pi — 1 
using (l2j 

™~l n c m— 1 n-m f n-2 \ 

^-("-^E ^: 1)!( :_ m) , (n(i-«o) f <» ^ 

m=l v ' K ' \j=m / 

Direct calculation gives leading and lowest powers of s in each of these polynomials. 

p n = n n ~ 2 s^-W"- 2 )/ 2 + ( n - 2) n"" 3 s ™(™- 3 )/ 2 + I ( n - 3) (7n - 6) n"" 4 sH"" 3 )" 1 

+ - (n - 1) (61n 2 - 338n + 384) n n ~ 5 S 5(*+i)(n-4) 
6 

+ 2~4 ^ ~~ ^ ( 705 ™ 3 ~~ 6265 ^ 2 + 17018/2 ~ 15 000) n n " 6 s |(™+ 1 )(«- 4 )- 1 + . . . 

+ Q (n - 1) (n - 2) 3 n ~ 3 - 1^ s 2 + ((n - 1) x 2 n " 2 - l) s + 1 . (36) 

So written, a power is understood to be absent when any exponent in its coefficient is < 0. 
There is also an exact "sum rule" for each polynomial. 

Pn (s = l) = ( (n-l)\) 2 . (37) 

The auxiliary function's leading asymptotic behavior is given for extreme values of s by 

oo n 

soE(^iv=^. (38) 

71=1 V 

* (*) s ^ £ („_!); = LambertW , (39) 

n=l ' ' 

as well as the formal result, 



*w a — v^ = tt —■ (40) 
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Recall in the power series for \l/ we chose to take the coefficient of x to be unity, but that 
normalization is arbitrary. So too are the normalizations for the asymptotic expressions. 

Now consider some particular functional roots and powers. An obvious check on the 
series for \l/ and is to verify that (s^> (x)) = sxe x , and indeed this is true to O (x 5 ) 
for the explicit results given in 021]) and (j25p . Also, when f\ (x) = sxe x it is well-known 
that the inverse function is /_i (x) = LambertW (-), as in ( I2T1) . We may check that the 
previous series for \l/ and do indeed give the series for LambertW when we take t = —1 
in f t (x) = (s* \I> (x)). For generic s, using the explicit series (jMJ) and (125]) . we find 

»- g» (,)) = i, + (-i) + ^ + (- A) x < + + o , (41) 

in perfect agreement with the textbook power series for LambertW (-x) ■ As another 
example, the functional square-root, such that f 1 / 2 (fi/2 (x)) = sxe x , follows immediately 
for any s. For instance, 

f 1/2 ( x ) = V' 1 (2tf (x)) = 2x + -x 2 + —x 3 + —x 4 ^—x 5 + O (x 6 ) , (42) 

ji/2k j v v J J ^ 3 45 135 137700 v ; ' v ; 

where we have avoided any irrational coefficients by the choice s — 4. This is quickly 
checked to satisfy fu 2 (fi/2 {x)) = 2xe x , to the order given. 

Integer iterates of fi (x)\ 8=1 = xe x are easy to construct, without approximation, and 
plot. We have fi (x) = xe x , f 2 (x) = xe x e X£X , etc. In general, we find the form 

fn+i (x) = ai exp (a x (1 + a 2 (1 + a 3 (• • • (1 + a n _i (1 + a n )) ■ ■ • )))) , (43) 

where a± = xe x , and for a n > 2 there is the recursion relation, 

a k+ i = exp ^II a ^ ' ^ 44 ) 

Plotting these reveals an ordered sequence of upward convex functions for x > —1, f n+ \ (x) > 
f n (x), and confirms that each has a minimum at fact easily established by the 

chain rule of differentiation. The minima are also ordered, f n+ i (— 1) > f n (— 1), with 
/i (—1) = — 1/e, and approach the x axis as n increases, lim^oo f n (—1) = 0. 

We plot four of these integer iterates. The curves are just time slices from a continuation 
of the ft (x) surface in Figure 10. All other positive iterates are also upward convex and 
ordered on the domain x > — 1, f t% (x) > f tl (x) for t 2 > t±, have f[ (— 1) = with minima 
at x = — 1 for t > 0, and are easily visualized as intercalated between the curves given in 
the following Figures. 
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Figure 11: First four integer iterates of xe x (fi, f 2 , fz, and / 4 ) plotted alternately in green 

and orange, above the identity map in black. 

Mirror imaging any one of these curves through the straight line of the identity map in 
the usual way gives the upper branch of the corresponding inverse function, e.g. /_i = 
LambertW, /_ 2 , f-3, and /_ 4 , and establishes the radius of convergence of the Taylor series 
expansion for /_ n (x) to be just R n = \ f n (— 1)|. 

We also plot some s — > 1 series approximations to the fractional iterates for — 1 < x < 0. 
These again give an ordered sequence of curves between fx and the identity map. 
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Figure 12: O (x 10 ) series approximations for fractional iterates fi (at top), /1/2, /1/4, /1/8, 
/1/I6, and /1/32, plotted alternately in orange and green, above the identity map in black. 

Once more these curves are just time slices from the f t (x) surface in Figure 11. For a 
measure of the accuracy of the O (x 10 ) series approximations used to make this last plot, 
note that f[ (—1) does not always vanish in the plot, but for the exact function, it should. 
In the plot, disparities are clearly discernible as x — > — 1, especially for the orange curves, 
where instead of vanishing we see that f[, 2 (— 1) > 0, f[/ 8 ( — 1) < 0, and fy 32 (~ 1) > 0- Nu- 
merically from the series, £f(x,l)\ x=l = -2.755 7 x KT 6 , £f (x, 1/2) = 0.19512, 
f-fOM/C^ = -0.03058, £7(^1/8)1^ = -0.29^32, fj £T/16) \ x= _ t = 
—0.030 586, and ^/ (x, 1/32) | __ = 0.346 87. These numerical disparities exist because 
x = — 1 lies precisely at the radius of convergence for the power series in question and is a 
branch point in the exact functions for generic t, so the series become poor approximations 
as x — > — 1. Better numerical results for fractional iterates near x = — 1 can be obtained 
by incorporating branch points into approximate trial solutions of the functional equation, 
and then matching these solutions onto the series for x > — 1. This is work in progress. 

V. SUMMARY FOR REPULSIVE POLYNOMIAL POTENTIALS 

The method of the paper can be used directly to solve for trajectories in well-known 
polynomial potentials, for motion towards or away from unstable fixed points, when turning 
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points are not encountered in finite time. Here we list the essential features for three 
examples. 



Quadratic: 



Quartic: [111 



Sextic: 



V (x) = -v 2 (x) = -x 2 (45) 

V (x) = X 

2^ (a:) = * (2a;) 
* (x) = x , (x) = X 

x (t) = xe l , t = r In v4 , x = 



V(x) = -v 2 {x) = -l + 2x 2 -x 4 (46) 
v (x) — (1 — x) (1 + a;) 



3 v ' \x + 2 

t /n fx — 1 \ T 1 / s 2 + a; 
^(x) = 2 1 , (x) - 



x + 1 J 2 — x 

x-l + (x + l)e 2t R 
X(t)= l-x+(x + l)e^ * = ^^ 



x 




a; 2 + 2x 4 - x b (47) 



a: 



r In a/2 , xo = 



To produce these examples, we considered right-moving, zero-energy configurations. For 
convenience we rescaled t, and then we constructed the exact series solution to Schroder's 
equation in the neighborhood of selected fixed points, xq- From the auxiliaries, 
we then recovered the continuous time iterates. That is to say, x (t) = f T (x — xo) = 
(s Tl J/ (x — xq)), with r oc t and xq as given above. 



VI. CONCLUSIONS 



In conclusion, we suggest taking a broader perspective and considering other points of 
view in dynamics that invoke Schroder's functional equation, s o ^ = ^ o /, or its inverse, 
the Poincare equation, x -1 os = /o As usual, s is just the simple multiplicative map, 
or change of scale, s : x — > sx, while / is a less trivial, but given function. For example: 
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i) Single trajectory maps: Invert a trajectory function x (t) to obtain the time t (x), at 
least for some interval in t, assuming x (0) = 0. Then consider Schroder's functional 
equation written as x( x l v ) — V X (t ( x ) ; v). Is there a cogent relation between the 
parameter v and the initial velocity? What is the physical meaning of the resulting 
function x? 

The answers are straightforward. Write the equation to place emphasis on the time depen- 
dence: x ( x [t) ; v) = vx(t;v). Then analyticity near t — requires x (0; v) = 0, and, if 

x' (0; v) 7^ 0, then v = . We also have x (t) = x l ( v X (*))> where additional v de- 

at t=o 

pendence of x, if any, is now implicit. So, dynamical evolution along a single trajectory is in 
this sense a functional similarity transformation acting on the initial velocity: x = x~ lovo X- 
As a consequence of any non-commutativity between x an d the simple multiplicative map, 
we will have x~ X ° v ° X 7^ v i an d the trajectory will evolve as a nonlinear function of t - 
with linear t dependence being just free particle behavior: x (t) = vt. We may again think 
of x as an auxiliary function defined on the initial phase space which encodes all solutions 
of the classical equations of motion. 

ii) Pseudo-scaling: Given a trajectory, again construct the inverse function x~ l . Then 
what is the significance of the new time variable T (t; s) = x~ l (sx (t))? 

This change of time variable simply rescales the solution. That is to say, X (T) = x (T (t)) = 
sx (t). Note that T (t) differs from a linear function of t only if the trajectory function fails 
to commute with the multiplicative map. 

iii) Iterative time evolution: From a lattice of time points, is it always possible to use 
Schroder's functional equation to obtain continuous time evolution? 

In this paper we have discussed how, under certain circumstances, the answer to this last 
question is affirmative. It remains to see whether the method can be applied in all situations. 
Nevertheless, we believe that this particular functional method will be applicable to problems 



in classical chaos [151 . Il6l |. to complement existing methods for analyzing unstable fixed 
points. We also think there is no difficulty, in principle, to prevent the method from being 
applied to classical dynamics in higher dimensions, or even to quantum systems, at least 
under certain circumstances, upon extension to more variables [17j. We look forward to a 
time when uses of functional evolution methods have become commonplace. 
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